Stability of vortex structures in quantum dots 
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We study the stability and structure of vortices emerging in two-dimensional quantum dots in high 
magnetic fields. Our results obtained with exact diagonalization and density-functional calculations 
show that vortex structures can be found in various confining potentials. In non-symmetric external 
potentials we find off-electron vortices that are localized giving rise to charge deficiency or holes in 
the electron density with rotating currents around them. We discuss the role of quantum fluctuations 
and show that vortex formation is observable in the energetics of the system. Our findings suggest 
that vortices can be used to characterize the solutions in high magnetic fields, giving insight into 
the underlying internal structure of the electronic wave function. 

PACS numbers: 73.21.La, 85.35.Be 



I. INTRODUCTION 

Vortices can occur in quantum systems which are set 
to rotate, for example by applying an external magnetic 
field or by mechanical rotationiS^ Using the Gross- 
Pitaevskii mean-field approach, Butts and Rokhsar 3 
found that in a gas of rotating bosonic atoms which are 
weakly interacting with a repulsive force between them, 
vortices may form in a crystal-like lattice, in much anal- 
ogy to patterns that emerge in rotating superfluid helium. 
These vortex solutions appear as holes (vortex lines) in 
the particle densities, where each single zero of the Gross- 
Pitaevskii wave function corresponds to a unit vortex. 
With increasing angular momentum, the bosonic cloud 
develops a fiat shape, with more and more vortices pen- 
etrating it. 

The analysis of the electronic structure of two- 
dimensional parabolic quantum dots 4 - 5 (QD) in high 
magnetic fields has recently shown that, surprisingly, vor- 
tices can appear also in fermion systems showing many 
similarities to the boson case££ High-field solutions of 
the density-functional theory revealed zeros in the elec- 
tron densities, with electron currents circulating around 
them. These were interpreted as vortex clusters^ Many- 
body techniques could uncover vortex formation at high 
magnetic fields 6 '- or at large rotation!, giving credence 
to this interpretation. With increasing magnetic field 
or rotation, successive transitions between stable vortex 
configurations were found. These vortices carry magnetic 
flux quanta similarly to the fractional quantum Hall ef- 
fect (FQHE), where the system can be approximated by 
Laughlin wave functions, which attach additional vortex 
zeros at each electron^ In quantum dots 6 - and quantum 
dot molecules^ vortices were found to be bound on the 
electron positions but, unlike in the Laughlin wave func- 
tions, additional vortices appeared between the electrons. 
A single vortex bound to an electron is a Pauli vortex, 
because it is mandated by the exclusion principle. In 
general, the number of vortices on top of each electron 
must be odd for fermions to have correct particle statis- 
tics. The additional vortices not bound to electrons give 



rise to rotating currents of charge and a charge deficiency 
down to zero electron density at the vortex centra. By 
going around a vortex the wave function gains a phase 
change of 2irl, where / is the winding number or vortex 
multiplicity. 

As characteristic for mean-field theories, the self- 
consistent solutions of the Kohn-Sham equations for 
fermions, as well as the bosonic equivalent, the Gross- 
Pitaevskii equations, break the symmetry of the quan- 
tum state. For Bose-Einstein condensates under rotation, 
the Gross-Pitaevskii mean-field results have been shown 
to emerge as the correct leading-order approximation to 
exact calculations^ For Coulomb-interacting fermions, 
the density-functional approach suffers from the rela- 
tively crude approximations for the exchange-correlation 
energy, as well as from the problems in using a single- 
configuration wave function^ It is thus important to 
compare the fermion mean field approach to a more rig- 
orous solution of the full many-body Hamiltonian. How- 
ever, this is a tedious enterprise: analytic solutions are 
out of reach for N > 3, and a numerically accurate direct 
diagonalization is also typically restricted to fairly small 
particle numbers. 

For azimuthal symmetry of the confinement, the parti- 
cle density retains this symmetry in the laboratory frame 
of reference. Therefore, to study the appearance of vor- 
tices in those systems, a rotating frame, 13 conditional 
wave functions, 6 or correlation functions have to be ex- 
amined. In contrast to the bosonic case^ the pair cor- 
relations are not very informative for fermions due to the 
disturbing influence of the exchange hole. Instead, the 
vortex solutions have been pinpointed cither in a pertur- 
bative approach^ or using a conditional wave function 
which fixes N — 1 electrons to their most probable posi- 
tions, and then calculates the wave function of the Nth 
electron? 6 * The latter method is not unproblematic either, 
since the vortices are not independent of the electron dy- 
namics. The vortex locations in the conditional wave 
function depend on the positions of the fixed electron co- 
ordinates as well as on the choice of the probing electron. 

In the FQHE the off-particle zeros are usually con- 
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trasted to on-particle zerosi 1 ^ The on-particle zeros are 
independent of electron coordinates except the coordi- 
nate of the electron to which the vortex is attached. On 
the other hand, the off-particle zeros of a given parti- 
cle are not necessarily off-particle zeros of other parti- 
cles. Graham et alw& used this fact to conclude that off- 
particle zeros are not vortices in the real sense and there- 
fore no charge deficiency is necessarily associated with 
them. Analogously to the FQHE case, the conditional 
wave functions indicate that the electrons in quantum 
dots see different positions for the off-particle vortices, 
which is a manifestation of the fluctuations in the sys- 
tem. 

In this paper we present direct theoretical evidence 
that in quantum dots only the off-electron zeros show vor- 
tex structures in the electron and current densities, and 
conclude that the vortex structures are stable and can 
be used to classify the internal structure of the many- 
electron wave function. We do this by applying non- 
symmetric confining potentials which leads to the local- 
ization of the vortices. In the rotationally symmetric 
potential vortices are not localized and they move as the 
electron coordinates are changed averaging out the effect 
of vortices on the particle and current densities. In ro- 
tationally non- symmetric potentials, however, the lower 
symmetry should cause (at least partial) localization of 
vortices directly in the particle density. The breaking 
of the circular symmetry was used already by Manni- 
nen et al. in the study of Wigner localization of elec- 
trons in elliptical QD's using exact diagonalization (ED) 
techniques jii Here we show that the same trick can be 
applied to vortices which become directly visible in the 
exact particle and current densities. We find, however, 
that the Pauli vortices at the electron positions do not 
contribute to this effect. In our calculations even a small 
asymmetry in the confining potential is sufficient to lo- 
calize vortices and the fluctuations do not destroy this 
effect. We compare the exact results to the mean field 
solutions and suggest that the results can be generalized 
to arbitrary geometries. 



II. SYSTEM CHARACTERISTICS 

We consider N electrons trapped by a confining po- 
tential V c and subject to a perpendicular, homogeneous 
magnetic field B = (0, 0, B). The system is described by 
an (effective-mass) Hamiltonian 



atomic units, i.e., in Ha* 



11.86 meV and in a* B ~ 9.79 



H = 
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+ V c (Ti) 



1 



4ne ; 

i<j 



(1) 

where A is the vector potential of the magnetic field B, 
m* the effective electron mass, and e is the dielectric 
constant of the medium. We apply the typical mate- 
rial parameters for GaAs, namely, m*/m e = 0.067 and 
e/eo = 12 A. We give the energies and lengths in effective 



At high magnetic fields, after complete polarization of 
the QD, the exchange energy results in the formation of 
a stable and compact structure, the so called maximum 
density droplet (MDD)>±& It is a finite-size equivalent of 
the v = 1 quantum Hall state which assigns one Pauli 
vortex at each electron position. The MDD state can 
be found in various QD geometries. In circularly con- 
fined QD's the electrons occupy successive angular mo- 
mentum states on the lowest Landau level, which in the 
parabolic case leads to a relatively flat electron density. 
In rotationally non-symmetric potential wells, instead, 
the MDD window can be deduced from the kinks in the 
chemical potentials. In addition, the magnetic field for 
the MDD formation can be accurately predicted from the 
number of flux quanta penetrating the QDiiS 

As the magnetic field is increased, the compact elec- 
tron droplet is squeezed and eventually the repulsive in- 
teractions between the fermions cause the MDD to re- 
construct. For parabolic QD's, different scenarios of 
the reconstruction have been suggested. Chamon and 
WerJS found a "stripe phase" where a lump of electrons 
separate s from the MDD at a distance rj 21b, where 
Ib = \J h/eB is the magnetic length. Goldmann and 
Renn introduced projected necklace states which they 
found to be lower in energy than the states found by 
Chamon and WeniSI Geometrically unrestricted Hartree- 
Fock 2 ^ and CSDFT 23 ' 24 studies suggested that such edge 
reconstruction would occur with a modulated charge den- 
sity wave along the edge. For a sufficiently small Zeeman 
gap, this polarized reconstruction may be preempted by 
edge spin texturesi2£i2& 

The ED shows instability with respect to addition of 
internal holes as discussed by Yang and MacDonald. 27 
These holes were recently re-interpreted as vortices in 
mean- field density- functional calculations^ They were 
found also further away from the dot center where the 
electron density is low. Vortices behave often like clas- 
sical (localized) particles in the mean-field approach and 
the vortices appear as rotating currents of charge with a 
zero in the particle density at the vortex centra. Several 
charge-density-wavc states that mix different eigenstates 
can also be interpreted as solutions describing a transport 
of a vortex to the center of a QD. 24 Formation of vortices 
causes usually a broken symmetry in the mean-field par- 
ticle density, even when the Hamiltonian is azimuthally 
symmetric. 

Tavernier et al.® studied distribution of zeros in the ex- 
act many-body wave function of systems containing up to 
4 electrons. They compared the results with the rotating- 
electron-molecule (REM) wave functions^ In the regime 
where the effect of the external confining potential can 
be neglected, the rotating electron molecule model can 
provide an intuitive description of the Wigner-localized 
electrons. Tavernier and coworkers^ found out, however, 
that the REM model is unable to predict the clustering 
of vortices near electrons. 
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III. NUMERICAL RESULTS 
A. Mean-field description of vortex solutions 

In order to solve the many-body Schrodinger equation 
corresponding to the Hamiltonian we first work in the 
mean-field picture and apply the spin-density-functional 
theory (SDFT). For the self-consistent solution of the 
Kohn-Sham equations we employ a real-space scheme^ 
where the external confining potential V c can be arbitrar- 
ily chosen without symmetry restrictions. The exchange- 
correlation effects are taken into account using the local 
spin-density approximation (LSDA)^ At high magnetic 
fields, the effect of currents in the exchange-correlation 
potentials becomes non-negligible, and the current-spin- 
density-functional theory (CSDFT) 31 gives a slightly bet- 
ter approximation to the ground state energy^ The CS- 
DFT is computationally more demanding than the SDFT 
but according to our test calculations qualitatively sim- 
ilar vortex structures were found to emerge in both for- 
malisms. We apply the SDFT throughout this paper, 
since we found in these tests that it captures all the es- 
sential physics of these systems at much lower computa- 
tional work. 

The confinement is chosen to be a two-dimensional har- 
monic oscillator potential with elliptic deformation, de- 
fined as 

V c (r) = ^ 2 (Sx 2 + ^y 2 ), (2) 

where x and y are the major axes of the ellipse, Hlo is the 
confinement strength, and 6 is the eccentricity. 

For comparison, we also a pply a rectangular hard-wall 
confinement, defined in Ref. [191 The parameter describ- 
ing the deformation of the confining potential in this case 
is the side- length ratio (3 — L x / L y . In the post-MDD do- 
main, the SDFT predicts the formation of vortices inside 
the QD's. This is visualized in Fig. ^ showing solutions 
containing up to three vortices in elliptic (a-d) and rect- 
angular (e-h) six-electron QD's. The eccentricity 6 and 
the side-length ratio /3 have been set to 2 in these systems, 
respectively. As the vortices repel each other, the MDD 
states shown in Figs, ^a) and (e) reconstruct into states 
that enclose a linear vortex pattern along the longest 
major axis. There is a remarkable qualitative similar- 
ity in the high magnetic field behavior of these systems. 
However, a linear vortex cluster requires a rather large 
eccentricity (side- length ratio) of the QD. In the ellip- 
tic case with a confinement strength of huj = 0.5 Ha*, 
the triple-vortex configuration changes from triangular 
to linear when S is increased to about 1.4. 

Figure^shows currents induced by the magnetic field. 
The current is flowing clockwise around the vortices and 
anticlockwise on the edges of the dot. The reversal of the 
current near the vortex core is due to inner circulation 
of the electrons^ The number of vortices in the QD 
increases with the magnetic field, and the current loops 
of the vortices start to overlap. This causes formation of 



giant current loops which comprise several vortices. 

Fig-Elshows the chemical potentials, /-i(A) = E(N) — 
E(N — 1), of a rectangular ((3 — 2) QD containing 
N = 7 ... 16 electrons. Interestingly, the regime beyond 
the MDD is characterized by periodic oscillations in /i as 
a function of the magnetic field B. The peak positions 
in the oscillations match with the transitions between 
adjacent vortex states and mark the emergence of addi- 
tional vortices one- by-one in the QD. This is presented 
for N = 6 in Figs, ^f)-(h). Note that the oscillations 
get stronger as the number of electrons increases. The 
origin of the oscillations lies in the large reduction of the 
Coulomb energy in connection with the vortex forma- 
tion and the coexistent pronounced localization of the 
electrons. The oscillations are also visible in the total 
magnetization, M — —dE to t/dB, as shown in Fig. 4 in 
Ref. [H The magnetization in rectangular and elliptic 
six-electron QD's are compared in Fig. |31 In elliptic 
QD's the oscillations are weaker and less regular than in 
their rectangular counterparts where the size of the dot 
is a constant. The soft- wall confinement makes the dot 
more flexible in minimizing the total energy. These re- 
sults are in accord with the exact diagonalization results 
by Goldmann and RennSi, who calculated the chemical 
potentials of QD's in parabolic and non-parabolic "coffee- 
cup" shaped confinements. The latter confinement has 
a hard wall similarly to our rectangular potential well. 
Goldmann and Renn found kinks in the chemical poten- 
tials and the kink sizes increase with the electron num- 
ber. Moreover, the kinks in the non-parabolically con- 
fined systems were found to be much larger than in the 
parabolically confined systems. 

The results above suggest that the vortex formation is 
a considerable energetic effect that could be detected in 
appropriate experiments. Oosterkamp et alSSL measured 
the Coulomb oscillations peaks (chemical potentials) in 
transport experiments for vertical QD's and observed ad- 
ditional phase transitions beyond the MDD. They found 
oscillations in chemical potentials and the amplitude of 
these oscillations increased with the electron number. 
These data are consistent with our calculations and could 
indicate vortex formation in quantum dots. A direct in- 
terpretation of their result is, however, difficult due to 
the unknown shape of the QD sample and its eventual 
sensitivity to disturbance in the experiment. Moreover, 
the experimental oscillations may also indicate other phe- 
nomena, such as the formation of a spin texture, for ex- 
ample 

Conditional wave functions can be introduced not only 
for the analysis of the exact many-body wave function, as 
described in Ref. 6, but they are also useful to analyze the 
SDFT results. We use an auxiliary single-determinant 
function of the Kohn-Sham orbitals which emulates the 
exact conditional wave function^ This allows a study of 
the SDFT solutions where the electron density may have 
several minima, but the vortices are not directly local- 
ized to fixed positions. This may be due to the mixing of 
several eigenstates as shown in Ref. |3 For instance, the 
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FIG. 1: Spin-density- functional- 
theory (SDFT) electron densities 
and currents in an elliptic (a-d) 
and rectangular (e-h) six-electron 
quantum dot in different magnetic 
fields. The potential parameters in 
effective a.u. are (Hlj, <5) = (0.5, 2) and 
(L x , L y ) — (2\/2iT, v2tt) in elliptic and 
rectangular dots, respectively. The 
increasing of the magnetic field leads 
to a formation of a vortex pattern 
beyond the maximum density droplet 
(MDD) solution (B = 12 T) in both 
geometries. 



electron densities for the SDFT states at 18 T and 21 T 
show no density zeros. However, the corresponding con- 
ditional single-determinant functions show vortices near 
the fixed electron ring (see Fig.^J). In this picture Pauli 
vortices can be seen on the fixed electrons, and additional 
vortices are found between the electrons. 

The electron densities of the elliptical 6 = 2 dot in 
Fig. 0] show localization of both electrons and vortices. 
The electron localization results in six density maxima 
which is consistent with earlier SDFT calculations of el- 
liptically confined QD's>i£ The additional vortices give 
rise to charge deficiency in the electron density as seen in 
the right panel of Fig. According to the SDFT results, 
the intensity of the vortex localization strongly depends 
on the magnetic field. In the above example, the localiza- 
tion of vortices is partial with pronounced minima in the 
electron density in solutions of up to two vortices. As 
the magnetic field is increased further, the localization 
becomes complete, i.e., the electron density vanishes at 
the vortex core. This reflects the decrease of the mag- 
netic length is. 

The results indicate that only the additional vortices 
between the electrons have a charge deficiency associated 
with them. Conditional wave functions 6 and the total 



electron densities calculated with the SDFT show that 
there are vortices also further away from the dot cen- 
ter. In this region the electron density is usually a tiny 
fraction of the maximum electron density. The quantum 
fluctuations smooth out the effect of these external vor- 
tices and they cannot be directly observed in the exact 
particle density, even though the rotational symmetry is 
broken. 



B. Exact diagonalization for an elliptic dot 

We compare now the above results obtained within 
SDFT to those of a direct numerical diagonalization of 
the many-body Hamiltonian matrix. We hereby focus on 
small particle numbers and magnetic fields high enough 
such that the description is to a large extent restricted to 
what in the isotropic case would correspond to the lowest 
Landau level. In order to display the internal structure 
of the many-body wave function, we break the spherical 
symmetry of the dot by applying the elliptical confining 
potential [Eq. We assume full polarization of the 

electron droplet, and neglect the Zeeman energy. 

For the deformed case, 8^1, where the total angu- 
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FIG. 2: Chemical potentials for iV-electron rectangular (/3 = 
2) quantum dots as a function of the magnetic field. The dot 
is expected to be fully spin-polarized. After the MDD window 
there are oscillations in the chemical potential which matches 
with the appearance of additional vortices (one-by-one) into 
the quantum dot. 
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FIG. 3: Magnetization of a rectangular (solid line) and ellip- 
tical (dashed line) six-electron quantum dot as a function of 
the magnetic field. The ellipse eccentricity 5 and the rectan- 
gle side-length ratio f3 are both set to 2. The numbers in the 
figure denote the correponding number of vortex holes in the 
electronic structure. 



lar momentum is not any more a good quantum number 
(while, however, we still have good parity), the most ap- 
propriate and efficient basis set spanning the Fock space 
is formed by the eigenstates to the single-particle part of 
the Hamiltonian Eq. (Q. These must be calculated nu- 
merically. We determine the M lowest ones by directly 
diagonalizing the single-particle part of the Hamiltonian 
in a basis consisting of a sufficient number of correspond- 
ing Fock-Darwin 34 states at 8 = 1, which are known an- 
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FIG. 4: SDFT solutions of six-electron, elliptically (8 = 2) 
confined quantum dots at different magnetic fields. The con- 
finement strength is set to hut = 0.5 Ha*. Left panel: Con- 
ditional single-determinant wave functions of the Kohn-Sham 
states. The fixed electrons are marked with crosses and the 
probing electron is the rightmost electron at the top. The 
contours show the logarithmic electron density of the probe 
electron and the grey-scale show the phase of the wave func- 
tion. The phase changes from n to — n at the lines where 
shadowing changes from the darkest grey to white. The vor- 
tices that cause charge deficiency in the center of the dot are 
marked with + signs. Right panel: Electron densities show 
vortex holes and partial localization of the electrons forming 
six density maxima. 



alytically. Once the single-particle basis is at disposal, 
the Fock states are generated by sampling over all possi- 
bilities to set N particles on these M. states. From this 
sampling, only those Fock states with defined parity and 
configuration energy^ less than a defined cut-off energy, 
are chosen for diagonalization. The cut-off energy was 
adjusted to restrict the number of Fock states (i.e. the 
matrix dimension) to be less than about 50000. We limit 
the single-particle basis dimension to M. < 44 and use 
numerical integration for calculating the Coulomb ma- 
trix elements. The many-body Hamiltonian Q is then 
diagonalized in the obtained subspace. Densities (n(r)) 
and real currents (j(r)) are finally calculated in order to 
compare the broken-symmetry solutions of the ED di- 
rectly to the mean-field results. Here, 



n(r) = y^5(r - rj 



(3) 



is the density operator. The real current is obtained by 
taking the expectation value of 



j(r) =j p (r) + — A(r)ra(r) , 



(4) 
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FIG. 5: Electron densities (grey scale) and current densities 
(arrows) of exact diagonalization (ED) solutions for an ellip- 
tical quantum dot with 8 = 1.1. (a) The ground state with a 
two- vortex structure, (b) The first excited state with a single 
vortex at the origin. The right column shows the correspond- 
ing electron densities at the longest major axis of the ellipse. 
The confinement strength is Hu) = 0.5 Ha* and the magnetic 
field is B = 13.4 T. 



where 



V i( 5(r - !•*)] (5) 



is the paramagnetic current operator. 

We should note at this point that for obtaining an ac- 
curate description of the total energy of the system, using 
only 44 lowest single-particle states is not sufficient for a 
full convergence of the total energy. However, the rela- 
tive energy differences and the geometrical structure of 
the electron and current densities of the ground-state and 
lowest-lying states were converged within this basis set. 
Since our aim is the comparison of the broken-symmetry 
many-body with the mean-field solutions, rather than a 
detailed discussion of energy spectra and excitation en- 
ergies, this truncation appeared reasonable. 

Figure [3] shows the electron and current densities, 
(n(r)) and (j(r)), of the ground state and first excited 
state for an elliptic dot with eccentricity 5 = 1.1 and 
magnetic field 13.4 T. Both states have parity ir = — 1, 
and the first excited state is separated from the ground 
state by only 7.7 mHa*. The ground state shows a vortex 
pattern around two minima in the density and its total 
angular momentum is L — —25.4ft. The leading single- 
particle configuration^ of the Fock state has the form 
| 1001111100...), with amplitude | c c | 2 = 0.4. The first 
excited state with L = —21 .4ft shows a pronounced sin- 
gle vortex at the dot center, with the current circulating 
clockwise around the origin. The hole shows the charac- 
teristic cone shape for a vortex and appears nearly local- 
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FIG. 6: Electron densities (grey scale) and current densities 
(arrows) for two- vortex solutions in elliptical QD's. (a) Ex- 
act diagonalization (ED) result at S — 1.1. (b) ED result at 
S = 1.2. (c) SDFT result at 5 = 1.2. The inset shows the cor- 
responding electron densities at the longest major axis of the 
ellipse. The vortices are localized in the SDFT solution, and 
the ED solution shows slight derealization due to quantum 
fluctuations. The confinement strength is hcu = 0.5 Ha* and 
the magnetic field is B — 17 T. 



ized at the dot center. The leading configuration of the 
Fock states is | 011111100...) with amplitude | c c | 2 = 0.7. 
The second excited state has parity ir = +1 and is sep- 
arated by a gap of 53 mHa* from the ground state. It 
has angular momentum L — —20.4ft and no clear vortex 
structure. 

Figure [S] shows the densities and density contours for 
the cases S = 1.1 (a) and S — 1.2 (b) at an increased 
magnetic field of 17 T. Both states have parity ir = — 1 
and angular momentum L = —26.7ft (a) and L = —28.5ft 
(b). In both cases, the current circulates around fairly 
pronounced minima in the charge density, representing 
two vortices at the dot center. When the eccentricity in- 
creases from 5 — 1.1 (a) to 6 = 1.2 (b), the electrons 
begin to arrange themselves in the form of a Wigner 
molecule with counterclockwise rotation of the current 
around the maxima of the electron density, and clock- 
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wise rotation around the vortices (i.e. the minima in the 
density) , as it was also observed at S = 1.1. The reversal 
of the current near the vortex core is consistent with the 
SDFT results. 

It is important to note here that the ground state in 
both cases is a state with tt = +1, which for (a) is 
34 mHa* and for (b) is 7.6 mHa* lower in energy. No 
states were found in between these two lowest states with 
7r = +1 and 7r — — 1, i.e. both states shown in (a) and 
(b) appeared as first excited states, respectively. The 
fact that the energetic sequence of the states changes 
with increasing field is, however, not surprising: It is 
well known 5 from the ED studies of circularly symmetric 
QD's, that the overall structure of the spectrum is inde- 
pendent of the magnetic field, the role of which is mainly 
to tilt the spectrum so that the minimum energy is at a 
different state (see Manninen et allL, and Maksym and 
Chakraborty 36 ). The above results from ED are com- 
pared with the mean-field calculation in Fig. (c) for 
5 = 1.2. Due to the approximations made in SDFT, the 
transitions between different vortex solutions appear at 
slightly different B-values. This fact, together with the 
tilt of the ED spectrum with increased B, can explain 
why in Fig. the first excited state in the ED solution 
compares to the ground-state two-fold vortex obtained 
by SDFT. In the ED result, quantum fluctuations de- 
stroy the complete localization of the vortices, i.e. the 
exact density at the vortex center still is about one third 
of the maximum density while in the mean-field result 
the density at the vortex center is reduced to zero. 

Increasing the eccentricity still further, the localization 
of electrons leads to the formation of a charge-density- 
wave-like crystal. Vortices in between the classical elec- 
tron positions become even more apparent. Figure [7| 
shows the ED electron densities (grey scale) and current 
densities (arrows) for larger eccentricity 8 = 2. The mag- 
netic field is (a) 11 T, (b) 17 T and (c) 22 T, and tt = -1 
in all cases. (The tt = +1 states are slightly higher in 
energy). The solution at 11 T shows the MDD, and vor- 
tices form at higher magnetic fields. These solutions can 
be compared to the SDFT solutions in Fig. ^ There is 
again a qualitative agreement, but now the ED shows 
stronger localization of electrons at both ends of the dot. 
While now the Fock states show much mixing and a clear 
dominance of a few single configurations could not be 
observed, the general trend is very similar to the SDFT 
result shown above in Fig. 



IV. SUMMARY 

We have found that in non-circular confinements the 
vortices can be directly seen in the exact diagonalization 
electron densities as holes or charge deficiencies around 
which the current is circulating. The vortices are sta- 
ble against the quantum fluctuations and even a slight 
asymmetry in the confining potential cause the vortices 
to show up in the exact many-body electron density as 
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FIG. 7: Electron densities (grey scale) and current densities 
(arrows) of ED solutions for elliptical quantum dots with S = 
2. The magnetic field is (a) 17 T, (b) 21 T, and (c) 22 T. 
The inset shows the corresponding electron densities at the 
longest major axis of the ellipse. The confinement strength is 
hu = 0.5 Ha*. 

density minima. However, due to quantum fluctuations, 
there are no zeros in the electron density. The spin- 
density-functional calculations are in accordance with 
these results. They predict analogous vortex formation 
also in rectangular hard-wall quantum dots suggesting 
that the results can be generalized to a wide variety of 
geometries. The chemical potentials of the rectangular 
dot show features in the energetics which suggest the 
possibility of a direct comparison to experiments. 

In the light of these results it is justified to speak of 
vortices as real quasi-particles. Vortices are not indepen- 
dent of the electron dynamics, but they can be used to 
characterize the solutions in high magnetic fields which 
gives insight of the underlying internal structure of the 
electronic wave function. Prediction of vortex formation 
in hard- wall potentials gives also credence to the assump- 
tion that the vortices are robust and largely independent 
of the chosen geometry of the system. Therefore they 
seem to be universal features of the physics of the two- 
dimensional interacting fermion systems in strong mag- 
netic fields above the maximum density droplet forma- 
tion. 

To conclude, we briefly compare these results to 
bosonic systems, where vortex formation in rotating 
Bose-Einstein condensates has been much discussed both 
theoretically and experimentally. There are apparent 
similarities between the bosonic and fermionic caset In 
the composite fermion model bosons can be turned into 
fermions (and the other way around) by attaching fic- 
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tious magnetic fluxes on top of the electrons. Vortex 
formation appears as a universal phenomenon of quasi- 
two-dimensional quantum systems suggesting a synthetic 
theoretical rationale behind the phenomenal 
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